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ABSTRACT 

A catalogue of 14453 radio-loud AGN with 1.4 GHz fluxes above 3.5 mJy in the red- 
shift range 0.4 < z < 0.8, has been constructed from the cross-correlation of the NVSS 
and FIRST radio surveys with the MegaZ-LRG catalogue of luminous red galaxies 
derived from Sloan Digital Sky Survey imaging data. The NVSS provides accurate 
flux measurements for extended sources, while the angular resolution of FIRST allows 
the host galaxy to be identified accurately. New techniques were developed for extend- 
ing the cross-correlation algorithm to FIRST detections that are below the nominal 
1 mJy S/N limit of the catalogued sources. The matching criteria were tested and 
refined using Monte-Carlo simulations, leading to an estimated reliability of ~98.3% 
and completeness level (for LRGS) of about 95% for our new catalogue. 

We present a new determination of the luminosity function of radio AGN at 
z 0.55 and compare this to the luminosity function of nearby {z ^ 0.1) radio 
sources from the SDSS main survey. The comoving number density of radio AGN 
with luminosities less than 10^^ W liz~^ increases by a factor ^ 1.5 between z = 0.1 
and z — 0.55. At higher lumiosities, this factor increases sharply, reaching values of 
more than 10 at radio luminosities larger than 10^^ W Hz"^. We then study how the 
relation between radio AGN and their host galaxies evolves with redshift. Our main 
conclusion is that the fraction of radio-loud AGN increases towards higher redshift in 
all massive galaxies, but the evolution is particularly strong for the lower mass galax- 
ies in our sample. These trends may be understood if there are two classes of radio 
galaxies (likely associated with the "radio" and "quasar mode" dichotomy) that have 
different fuelling/triggering mechanisms and hence evolve in different ways. 

Key words: galaxies: evolution - galaxies: active - radio continuum: galaxies. 



1 INTRODUCTION 



Ever since the first studies of iLongaiij ||1966D . it has been 
widely accepted that evolution is required to explain the ra- 
dio source counts. Strong cosmological evolution of the high 
luminosity radio source population has been estab lished out 
to at least 2 = 2 — 3 fe.g. lDunlop fc Peacockiri990l ). with the 
space-densities of the most powerful radio sources increasing 
by a factor ~ 10^ relative to their densities at z ~ 0. How- 
ever, the evolution of the low-lum inosity radio galaxy pop - 
ulation remained more uncertain. IClewlev fc Jarvi i (l2004l ) 
found that the low luminosity population of radio sources 
with I/325MHZ < 10^^ W Hz~^sr~^ exhibits no evolution out 
to 2 ~ 0.8. In contrast. iBrown et all (|200ll ) found significant 
cosmic evolution out to z ~ 0.55, which they modelled using 
a pure luminosity evolution model as P{z) — P(0)(1 + z)*^ , 
with 3 < fci, < 5. Mo r e rece ntly, similar results were ob- 
tained by ISadler et al.l (|2007l ). who analyzed radio galax- 
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ies in the '2dF-S DSS Luminous Red Galaxy and Quasar' 
(2SLAQ) survey (ICannon et al.ll2006l ). They found that the 
number density of radio sources below Pi.4ghz 
Hz~^ grows by a factor of ~ 2 between z = and z = 0.55. 
These results confirm that the cosmic evolution of low power 
radio sources is considerably weaker than that of their higher 
radio power counterparts. However, because of the relatively 
small sizes of samples, none of these studies was able to accu- 
rately constrain the luminosity dependence of the evolution. 

Over the same period, our physical understanding of ra- 
dio galaxies has developed steadily. Fanaroff & Riley (197^) 
showed that radio galaxies exhibit a dichotomy in radio mor- 
phology. Low radio luminosity sources are usually 'edge- 
darkened', with radio jets that quickly decelerate and flare 
as they advance through the interstellar medium (Fanaroff- 
Riley class I sources, or FRIs). In contrast, most higher lu- 
minosity sources have jets that remain relativistic over kpc 
and up to Mpc scales, ending in bright hotspots ('edge- 
brightened' sources, or FRIIs). The dividing line between 
FRIs and FRIIs occurs around a monochromatic power 
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of log io(Pi.4GHz)~ 25 W Hz^^, although iLedlow fc Owen! 
l|l996l ) showed that it occurs at higher radio luminosity in 
more optically-luminous galaxies. This division is close to 
the break in the radio luminosity function, and has led 
to the hypothesis these two populations of radio galaxies 
correspond to the different populations of cosmologically- 
evolving sources, FRIIs being strongly evolving and FRIs 
much less so (e.g. Ijackson fc Walllll999l ). 

A more important division in the radio source popula- 
tion, however, appears to relate to the radiative properties 
of the AGN. Powerful F RIIs and a small subse t of the more 
powerful FRIs (e.g. see iHevwood et "al] l2007t ) show strong 
emissi on lines and e ither a visible or hidden quasar-like nu- 
cleus ijBarthellllOSgl ). indicative of radiatively efficient accre- 
tion. In contrast, most FRIs and a sub-population of weak- 
emission-lined low luminosity FRIIs show no evidence for a 
powerfully radiating nucleus (see Hardcastlc ct al. 2007 and 
references therein) suggesting radiatively inefficient accre- 
tion. This implies that a fundamental difference exists be- 
tween high and low power radio sources that is distinct from 
the radio morpological dichotomy. Instead, it is likely to be 
related to different origins or mechanisms of gas fuelling. 

Understanding this low luminosity radio source popu- 
lation, and its cosmic evolution, is not only important for 
understanding the physical origin of the radio activity, but 
also for understanding galaxy evolution. There is growing 
evidence that low-luminosity radio-loud AGN play a crit- 
ical role in regulating the masses and star formation rates 
of galaxies. Semi-analytic rnodels of galaxy forraation (for 
example see lWhite fc Frenk|[l99ll and lCole et"al]|2000l ) suc- 
cessfully reproduce many of the the observed properties of 
galaxies. However, these models over-predict the abundance 
of galaxies at the bright end of the luminosity function. This 
problem arises because many massive galaxies are predicted 
to sit at the centre of hot hydrostatic haloes, in which cool- 
ing flows are expected to develop. If the gas cools and forms 
stars, central group and cluster galaxies are too massive by 
the present day and hav e much bluer colours than observed. 
iTabor fc BinnevI (|l993l ) first suggested that radio galaxies 
could in principle regulate these cooling flows, preventing 
significant accretion of gas and limiting the mass of galax- 
ies. This concept of AGN feedback has now bee n successfully 
incorp orat ed into the semi-ana lytic models of iBower et al.l 
l|200gi and lCroton et al.l (|2006l ). 

Observational support for a scenario in which radio 
AGN play an important role in regulating the growth of 
massive galaxies has slowly been accumulating. Both lu- 
minous FRII and low power FRI sources can in principle 
inject a significant amount of energy into the surrounding 
gas, as their radio lobes expand and interact with the sur- 
rounding medium (~ 10^^ ergs~^ fo r the jets of a borderline 
FRI/FRII object, see[Bicknd][l99i). This has been directly 
observed in clusters of galaxies, where the expanding radio 
lobes create buoyan t bubbles and drill cavities in the hot in- 
tracluster medium (iBohringer et al.] 1 19931 : iMcNamara et al] 
I2OO0I : iFabian et aLll2003l ). In nearbv systems, which can be 
studied in detail, the energy supplied by the radio source is 
well-matched to t hat required to balance the cooling (e.g. 
iFabian et al]|2003l ). 

With the availability of large redshifts surveys like 
the 2-degree-Fiel d Galaxy Redshift Survey (2dFGRS, 
IColless et aLll200ll ) and the Sloan Digital Sky Survey (SDSS, 



lYork et aLlbood ) it has become possible to obtain the sky 
and redshift coverage needed to define sufficiently large sam- 
ples of nearby radio sources to study population statistics 
and global energetics. iBest et al] (j2005a) constructed a sam- 
ple of 2215 radio-loud AGN brighter than 5 mjy with red- 
shifts 0.03 < z < 0.3 using the SDSS Data Release 2 (DR2) 
in combination with the National Radio Astr onomy Qbser- 
vatory Very Large Array Sky Survey (NVSS. [Condon et al.l 
1 19981 and the Faint Images of the R a dio Sky at Twenty 
centimeters (FIRST. iBecker et al.|[l995l ) . iBest et"al] l|2005bh 
found that the fraction of radio-loud AGN, fnL, is a strong 
function of black hole and stellar mass, with dependencies 
fnL oc M^'^ and Jul oc Mgfj. Twenty-five percent of the 
very most massive galaxies host radio sources more power- 
ful than 10'^'' W Hz~^. By combining these resul ts with es- 
timates of the energetic output of these sources, iBest et all 
(2006) showed that the time-averaged mechanical energy 
output of radio sources in elliptical galaxies is comparable 
to the radiative cooling losses from the hot halo gas, for el- 
lipticals of all masses. Radio sources can therefore provide 
sufficient energy to suppress cooling flows and regulate the 
growth of massive galaxies in the nearby Universe. 

The role of powerful (radiatively-eflicient) FRII ob- 
jects could not be addressed by Best et al. due to the 
scarcity of such systems in the low redshift catalogue. An- 
other open question is how the heating-cooling balance in 
massive galaxies evolves with redshift. In order to address 
these issues, we require a large radio galaxy sample at higher 
redshifts. In this paper we extend the work of Best et al. 
by constructing a comparable catalogue of radio-loud AGN 
over the redshift range 0.4 < z < 0.8. The catalogue con- 
tains a large number of luminous radio sources. By compar- 
ing with a sample of local radio sources selected from SDSS 
DR4, we aim to study not only of the evolution of the radio 
galaxies, but also the evolution of the duty cycle of radio 
AGN activity in galaxies as a function of their stellar mass. 

This paper is organized as follows. In Section 2 we will 
describe the surveys and samples used in this work, as well as 
the matching procedure. In Section 3 we will present results 
on the evolution of the radio luminosity function, radio-loud 
AGN fractions and the bivariate radio luminosity-stellar 
mass function, and the stellar mass and colour distributions 
of radio AGN host galaxies. Finally, we will summarize our 
results in Section 4 and discuss the implications of our work. 

2 THE RADIO-OPTICAL GALAXY SAMPLES 

2.1 The MegaZ-LRG Galaxy Catalogue 

The Sloan D i gital Sky Survey l|York et all I2OO0I : 
IStoughton et al] |2002| ) is a five-band photometric and 
spectroscopic survey that has mapped almost a quarter 
of the whole sky, providing precise photometry for more 
than 200 million objects and accurate redshifts for about 
a millio n galaxies and quasa rs. Imaging data in the ugriz 
bands (|Fukugita et alj 1 19961) has been obt ained with a 
large format CCD camera l|Gunn et al.lll998l ) mounted on 
a dedicated telescope located at Apache Point Observatory 
in New Mexico. Astrometric coordinates are accurate to 
~ 0.1 arcsec per coordinate. Throughout this work, unless 
otherwise stated, all magnitudes are model magnitudes. 
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derived from the best-fit de Vaucouleurs or exponential 
profile in the r band. The amplitude is then scaled to 
fit measurements in other filters. This provides the best 
estimate of galaxy colours since the same aperture is used 
in all passbands. A correction for foreground Galactic 
extinction is applied to a l l the magnitudes following 

ISchlegel. Finkbeiner fc David (Il998l). 

The MegaZ-LRG feo Uister et al]|2007f ) is a photomet- 
ric redshift catalogue based on imaging data from the fourth 
Data Release (DR4) of the SDSS. It consists of ~1.2 mil- 
lion Luminous Red Galaxies (LRGs) with limiting magni- 
tude I < 20 over the redshift range 0.4 < 2 < 0.8. MegaZ 
adopts various colour and magnitude cuts to isolate red 
galazies at 0.4 < z < 0.8. The cuts are very similar to 
those ad opted by the '2dF-S DSS LRG and Quasar' project 
f2SLAQ. ICannon et al.ll2006l ). The specific colour cuts are: 



0.5 < 5 - r < 3 

r — i < 2 

cpar = 0.7{g -r) + 1.2(r - i - 0.18) > 1.6 



{r-i) 



r)/8 > 0.5 



(1) 
(2) 
(3) 
(4) 



Essentially, Cpar isolates early-type galaxies (LRGs), 
while dperp selects galaxies with redshifts z > 0.4. Accu- 
rate photometric redshifts are available for the entire LRG 
sample. These are derived using a neural network photo - 
metric redshift estimator fANNz. ICoUister fc Lahav|[20o3 ). 
The neural net was trained using a sample of ~ 13000 LRGs 
with spectroscopic redshifts from 2SLAQ. The r.m.s. average 
photometric redshift error for all the galaxies in the sample 
is arms = 0.049. 



2.2 The NVSS and FIRST Radio Catalogues 

The NRAO VLA Sky Survey (NVSS: ICondon et al.|[l998l ) 
is a radio continuum survey at 1.4 GHz, covering the whole 
sky at declinations (5 ,72000 > —40° . The survey provides total 
intensity and polarization measurements down to a limiting 
point source flux density of ~ 2.5 mjy. The rms noise and 
confusion limit in the total intensity maps is oaMS ~ 0.45 
mJy beam~^, except near very strong sources. The survey 
was carried out using the VLA in compact D and DnC con- 
figurations, which provides an angular resolution of around 
45 arcsec. A catalogue of ~ 1.8 million point sources was de- 
rived by fitting elliptical gaussians to the discrete peaks in 
the images. The catalogue lists peak flux densities and best- 
fit sizes for all sources with fluxes Speak ^ 2 mJy beam~^. 
The completeness of the NVSS survey is around 50% at 
~ 2.5 mJy, but rises sharply to more than 99% at 3.5 mJy. 

The VLA Fai nt Images of the Ra dio Sky at Twenty Cen- 
timeters (FIRST; iBecker et al]|l995l ') is a survey that covers 
10000 deg'^ of the North Galactic Cap, using the VLA in its 
B configuration at 1.4 GHz. This delivers a resolution of ~5 
arcsec and a typical noise of a rms ~ 0.13 mJy, except in 
the vicinity of bright sources (> 100 mJy) where sidelobes 
can lead to an increased noise level. The sky coverage over- 
laps most of the SDSS survey in the Northern Galactic Cap. 
As with the NVSS, a source catalogue including peak and 
integrated flux densities and sizes was derived from fitting 
a two-dimensional Gaussian to each source generated from 
the co-added images. Astrometric accuracy in this source 



catalogue is better than 1 arcsec for the faintest sources de- 
tected at the nominal threshold of 1 mJy. 



2.3 Identification of radio— loud AGN 

The NVSS and FIRST surveys are highly complementary in 
many respects. NVSS has sufficient surface brightness sen- 
sitivity to provide accurate flux measurements of extended 
radio sources with lobes and jets. On the other hand, the su- 
perior angular resolution of FIRST is necessary to identify 
the central core component of the radio sources and asso- 
ciat e it with th e host g alaxy detected in optical surveys. 

iBest et al] (|2005a ) developed an automated hybrid 
cross-matching method to identify radio galaxies in the 
SDSS DR2. The method made use of both NVSS and 
FIRST, thereby exploiting the advantages of each survey 
and avoiding the errors that arise by using them sepa- 
rately. A collapsing algorithm was used to identify multiple- 
component FIRST and NVSS sources. The final catalogue 
consisted of 2712 radio galaxies brighter than 5 mJy (2215 
radio-loud AGN and 447 star forming galaxies). The relia- 
bility and completeness of the catalogue was estimated to 
be 98.9% and 95%, respectively. 

Th is paper presents an extension of the iBest et al] 
(|2005al ') method by combining NVSS and FIRST with the 
high redshift MegaZ-L RG catalogue. We largely follow the 
procedure described in lBest et al.l (|2005ah , but a few small 
modifications were implemented in order to accommodate 
the characteristics of the new sample. The changes were 
tested and progressively refined using a set of 10 Monte- 
Carlo simulations of a 720 deg^ patch of the survey area. 
The numbers quoted in the description below pertain to this 
small test area. 



(i) iBest et al. 1 l|2005ah considered radio sources above a 
limit of 5 mJy. The major reason for this choice was to 
reduce any loss of sensitivity to multi-component NVSS 
sources. The MegaZ radio sources are more distant and their 
angular sizes are thus smaller, so there are fewer multi- 
component sources (see Section 2.5). In order to probe 
fainter radio galaxies out to higher redshifts, we adopted 
a fiux limit of 3.5 mJy. This limit corresponds to approxi- 
mately 7 times the noise level in the NVSS maps. 

(ii) Figure 1 shows the distribution of positional offsets 
between the SDSS LRG and the nearest NVSS match, for 
galaxies not classified as multi-component NVSS sources. 
There is a clear excess of galaxies compared to random 
matches for separations less ^10-15 arcsec. The excess re- 
mains significant out to ^^70 arcsec. A significant part of 
the effect at larger angular separations is bec ause galax- 
ies are intrinsically clustered (see discussion in iBest et al.l 
(|2005al ) and also Section 2.5). However, some part may 
consist of true associations between LRGs and extended 
NVSS sources. The subseq uent FIRST cross-r natching helps 
us pinpoint such cases. iBest et al] (|2005al ) carried out 
this cross-matching for all SDSS galaxies with a single- 
component NVSS match within 30 arcsec. In this work, this 
limit is increased to 60 arcsec. This identifies an extra 37 
(2.8%) genuine radio galaxies at a cost of 0.8 false matches, 
i.e. a decrease in reliability of only ~0.01%. 

(iii) Because the N VSS radio s o urces r each a fainter limit 
than those studied bv lBest et all (|2005al ). we find more ob- 
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Figure 1. Distribution of offsets between the MegaZ-LRG galax- 
ies and the closest NVSS source (solid line) and the corresponding 
offset distribution derived from the random samples (dot— dashed 
line). Most true matches clearly stand out at offsets less than 15 
arcsec, but the excess remains significant out to separations of 
~70 arcsec. 

jects (~8%) with no FIRST match. These might either be 
extended sources missed by the low surface brightness sen- 
sitivity of FIRST, or because the radio source has faded 
between observations. A careful inspection of FIRST im- 
ages (along with NVSS and SDSS cutouts) around the posi- 
tions of these objects revealed that in many cases there were 
bright spots near or slightly below the survey flux limit, 
that were clearly associated with a MegaZ optical galaxy. 
A method to use this information to increase the number of 
matches was developed and tested. The details are explained 
in the following subsection. 



2.4 Improved matching of sources without FIRST 
counterparts 

The high resolution of the FIRST survey is essential to al- 
low reliable cross-matching with dense optical samples. The 
NVSS sources without FIRST counterparts are the largest 
source of incompleteness in the final sample (see lBest et al.l 
l|2005al )). However, blindly accepting all LRG candidates out 
to a separation of 20 arcsec from the NVSS source introduces 
an average of 42.2 false detections in our random samples, 
reducing the reliability of the catalogue to 79% within our 
test area (see Figure 2, where the distribution of NVSS- 
MegaZ offsets are plotted for real and random data). 

We turned to the analysis of the FIRST radio maps 
around the candidate positions, in order to retain as many 
of these sources as possible without introducing significant 
contamination. Figure 3 shows FIRST cutouts of typical ob- 
jects with no catalogued FIRST source, but with a nearby 
NVSS-MegaZ association. In many cases one or more bright 
"spots" matching the NVSS source and MegaZ galaxy can 
be seen, while in some cases two well-defined symmetric ra- 
dio lobes centered around the optical galaxy are evident. 

We took advantage of this additional information in de- 
ciding whether candidates should be accepted as true FIRST 
matches. For galaxies within 10 arc sec of the NVSS so urce 
we accepted all matches, following iBest et all (|2005al ). In 
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Figure 2. Distribution of offsets between MegaZ— LRG galax- 
ies and NVSS candidates for sources that have no match in the 
FIRST catalogue. Results are shown for real (solid line) and ran- 
dom (dot-dashed line) samples. Matches with separations larger 
than ~15 arcsec are mostly the result of galaxy clustering. 

the test area, this provided 96 associations of which 9.9 are 
expected to be false (i.e. 10.3% contamination). For galaxies 
with offsets between 10 and 15 arcsec, we accepted matches 
according to the following procedure: 

(i) The FIRST image was segmented using an initial 
threshold of 3(Trms(a,5), where (Trms{ot,5) is the local rms 
noise in the map. This provided a number of seed spots 
whose centroids were refined using the intensity-weighted 
moments in both directions. 

(ii) A brightness profile was derived for each spot and 
both the half light radius and the flux inside this radius 
were determined. 

(iii) Any bright spot with a mean flux density greater 
than 0.4 mjy and within 5 arcsec of the optical counter- 
part was considered to be a valid match. These sources were 
accepted unless there were more than 4 bright spots inside 
a circle of 30 arcsec, in which case the match was rejected 
because it is then impossible to associate the source with 
a single MegaZ galaxy. Visual inspection of SDSS images 
confirmed this. 

To test this 'spot matching' technique we repeated the 
same procedure for matches without FIRST components in 
the random samples. This constitutes a robust check, since 
we are in effect looking at random locations in the sky, so 
we can assess how probable it is to find a spot pattern by 
chance. In the sky test area, our procedure identified 13 
additional FIRST for sources with NVSS-MegaZ offsets be- 
tween 10 and 15 arcsec, at the cost of introducing only 2.4 
false detections. This increases the overall completeness of 
the sample, for little cost in reliability. Our tests indicate 
that it is not possible to extract more matches beyond a 
~15 arcsec offset with reasonable confidence. 

2.5 Reliability, Completeness and Final Catalogue 

Table 1 lists the number of accepted radio matches for dif- 
ferent classes of radio sources found in the sky test area. We 
provide results for both the real and the random samples. 
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Figure 3. FIRST cutouts around sources without FIRST counterparts. The dashed circle marks the NVSS source position and size, 
while the plus sign indicates the MegaZ optical galaxy location. The small circles mark the detected spots. Cases (a),(b) and (c) are 
examples of accepted matches and cases (d),(e) and (f) are examples of rejected ones. A boxcar sliding average filter has been applied to 
enhance the visibility. 



as well as reliability estimates. Note that even though our 
sample is much deepe r and m ore than five times larger than 
the one of lBest et al.l (|2005al ). we have been able to compile 
a catalogue with comparable reliability (~98.3%). 

Because the true number of radio sources is unknown, 
the completeness is more difficult to estimate than reliabil- 
ity. Nevertheless, a reasonable estimate of the completene- 
ness can be obtained by considering the MegaZ-LRGs that 
have an NVSS source within a distance of 15 arcsec. These 
can be assumed to be a mixture of true matches, random 
alignments from sources at different redshifts, and a contri- 
bution from sources that are associated with other galaxies 
which just happen, through clustering, to be close to the tar- 
get galaxy. By using the Monte-Carlo simulations in the sky 
test area to estimate the latter two contributions, we can 
estimate the number of true matches and hence the com- 
pleteness. 

In the sky test area there are 1765 MegaZ-LRGs with 
an NVSS source within a distance of 15 arcsec. In compar- 
ison, there are an average of 225 sources within the same 
distance to the galaxies in the random catalogues. Deter- 
mining the contribution from clustered galaxies galaxies is 
more complicated. We attempted this in two ways. 

First, we selected all NVSS-SDDS candidates with a 
single FIRST component within 15 arsec, but which were 
rejected by the matching algorithm. We then searched for 
other optical MegaZ galaxies in a 3 arcsec region around 
these FIRST radio positions (a standard search radius for 
FIRST matching, e.g. .Obric et al. 2006). From a total of 
359 LRGs, 40 (11.1%) had a nearby LRG with an associ- 
ated FIRST source. The corresponding fraction in the ran- 
dom sample was 0%. This means 11.1% of the 1765 sources 
(~200) within 15 arcsec are due to clustering. 



To check this result, we estimated the excess number 
of LRG-LRG pairs over random b y integrating the angular 
correlation function w{6) given bv lRoss et al.l j2003) for the 
2SLAQ catalog out to 15 arcsec. This yielded value of 1.66, 
i.e. ~ 150 extra sources due to clustering. These two values 
are in rough agreement; the fact that the first method gives 
a slightly larger number could be explained by the fact that 
radio galaxies are more massive than the underlying LRG 
population, so we would expect them to be more strongly 
clustered. For our final completeness estimate, we will as- 
sume that 180 sources are due to clustering. 

The expected number of true matches is then esti- 
mated as 1765 — 225 — 180 = 1360. Our matching proce- 
dure detects 1309 objects, of which 21.7 are predicted to 
be false detections. The completeness can therefore be esti- 
mated as 100 * {detected — false) / expected or 100 * (1309 — 
21.7)/136 ~ 95%. This value is comparable to that ob- 
tained bv lBest et al.l (|200 5a). for a much lower redshift sam- 
ple. We note that this completeness estimate applies only 
to the MegaZ-LRG sample and does not take into account 
missing radio-loud QSOs or radio galaxies that do not meet 
the MegaZ-LRG colour cuts. Q 

Table 2 lists the final number of radio galaxies found 
in the MegaZ-LRG sample and included in the catalogue. 
Our catalogue is by far the largest sample of radio galaxies 
of its kind. Out of a total of 14553 sources, the majority 

^ A preliminary comparison with a complete spectroscopic cata- 
logue at similar redshifts indicates that ^ 20% of all galaxies with 
Af* > 10^^ M0 are too blue to fall within the MegaZ-LRG colour 
coundaries. Since radio AGN are biased towards the higher mass 
galaxies, the corresponding percentage should be smaller for this 
population. 
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Table 1. Numbers of different classes of accepted MegaZ-LRG 
radio galaxies, as well as the corresponding numbers found in the 
random samples. The numbers are for a 720 deg^ patch of the 
total survey area. Values in parentheses indicate the number of 
sources accepted using the spot-matching technique (see text). 



NVSS type 


FIRST type 


MegaZ 


Random 


Reliab. % 


NVSS sing. 


comp. 


109(13) 


12.3(2.4) 


88.7 


NVSS sing. 


1 comp. 


1068 


7.8 


99.3 


NVSS sing. 


2 comp. 


88 


1.2 


98.7 


NVSS sing. 


3 comp. 


38 


0.3 


99.2 


NVSS sing. 


4+ comp. 


6 


0.1 


98.3 


NVSS mult. 




24 


0.8 


96.7 


Total 




1333 


22.5 


98.3 



(78.6%) have single counterparts in both FIRST and NVSS. 
A significant fraction (11.5%) are resolved by FIRST into 
multiple components. Only 253 (1.7%) sources have multiple 
NVSS counterpart s, to be compared with the 6.2% found by 
iBest et al. 1 (12005^ 1. This is again expected, because the radio 
galaxies lie at larger dista nces. 

We note that as for iBest et all l|2005al ). a fraction of 
the candidate sources (1093 objects, or 7.5% of the sam- 
ple) with multiple NVSS or FIRST components were flagged 
as being too unreliable to classify by automated means. 
These were classified visually by simultaneous examination 
of SDSS, NVSS and FIRST images. Of these, 221 were ac- 
cepted as genuine sources. Finally, to ensure that the same 
NVSS source was not associated with two different MegaZ- 
LRG galaxies, the complete catalogue was checked. Such 
duplications arose for 560 candidates and all cases were vi- 
sually examined. Altogether, 564 positive matches (3.9%) 
result from the visual classification process. 

Finally, we note that it is important to distinguish be- 
tween stars and galaxies when compiling the radio galaxy 
catalogue, because late M-type stars have colours t h at are 
very similar to those of faint LRGs. ICollister et al] ()2007l ) 
used the same neural network technique used to derive pho- 
tometric redshifts, to perform star/galaxy separation. The 
neural net was given 15 different photometric parameters 
and trained using the 2SLAQ spectroscopic sample. The 
MegaZ catalogue provides a parameter, 5sg, which is the 
probability that a particular object is a galaxy rather than 
a star. In the subsequent analysis, we only consider objects 
with 5ag > 0.7 (99.93% of all the matched radio galaxies, 
96.78% of all optical galaxies). 



2.6 Comparison Sample 

In the following sections, we will compare the results ob- 
tained for the MegaZ-LRG catalogue of radio AGN to 
those obtained for galaxies selected from the main spec- 
troscopic galaxy sample, which have redshifts in the range 
0.01 <z < 0.30, with a mean redshift of {z) ~ 0.14. This is 
an extension to the fourth data release of the original DR2 
radio-matched catalogue of iBest et al] l|2005al ). 



Table 2. Total number of MegaZ-LRG radio galaxies with 1.4 
GHz flux above 3.5 mjy in the different classes. 



NVSS type 


FIRST type 


MegaZ 


NVSS sing. 


comp. 


1176 


NVSS sing. 


1 comp. 


11445 


NVSS sing. 


>1 comp. 


1679 


NVSS mult. 




253 


Total 




14453 



3 RADIO-LOUD AGN EVOLUTION 

3.1 Masses and (if + e)— corrections 

To derive stellar masses for each galaxy we used the kcorrect 
algorithm (Blanton fc Rowei s 2007), which finds the non- 
negative linear combination of spectral templates that best 
matches the fiux measurements of each g alaxy in a sense. 
These templates are based on a set of iBruzual fc Chariot! 
(2003') models and span a wide range of star formation his- 
tories, metallicities and dust extinction, so they can be used 
to estimate the mass-to-light ratio of a galaxy. This method 
yields stellar masses that differ by less than 0.1 dex on av- 
erage from estimates using other techniques, for example 
using the 4000A break strengt h and HS absorption in dex to 
estimate M/L as proposed bv lKauffmann et al.ll2003l . 

In order to compare radio galaxies at different redshifts 
in a fair way, we need to compute the volume over which the 
MegaZ sources are visible, given the magnitude and colour 
cuts of the surveys as well as the 3.5 mJy radio fiux limit. 
In this work we generated (K + e)-corrections using a li- 
brary of lBruzual fc CharlotI l|2003l) evolutionary stellar pop- 
ulation synthesis models. LRGs constitute a nearly homoge- 
neous population of old, red, early-type galaxies with spec- 
tral energy distributions that are reasonably well described 
by a single stellar population (SSP) model that forms all its 
stars in an instantaneous burst at hi gh redshifts (we ch oose 
z = 9.84) and then passively evolves (Wake et al.ll2006l ). We 
use these models to compute if -corrections, which account 
for the redshifting of the spectra through the different rest- 
frame bandpasses, and e-corrections, which account for the 
intrisic evolution in the spectral energy distribution of galax- 
ies because stars are younger at higher redshifts. With these 
corrections in hand, we are able to predict the u,g,r,i,z 
magnitudes as well as the radio luminosities of each MegaZ 
LRG as function of redshift, which then allows us to de- 
termine the precise minimum (zmin) and maximum (zmax) 
redshift that a particular object enters and then leaves the 
radio-optical sample. The weight for each galaxy is then 
l/Vmax, where Vmax is the integral of the comoving volume 
between Zmin and Zmax- 

As a check, we have computed the stellar mass function 
of the LRGs in the MegaZ catalogue and we compare our 
results with stellar mass function of low redshift massive 
galaxies drawn from DR4. This is shown in Figure 4. The 
MegaZ-LRG sample probes the high mass end of the stellar 
mass function, spanning a range in stellar mass from 10^^ 
to 10^^' ^Mq. Over this range in mass, there is essentially no 
evolution in the stellar mass function out to z 0.55. This 
is in accor d with recent res ults presented by iBrown et al.l 
l|2007l ) and I Cool et al.l l|2008h . who show that the upper end 
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Figure 4. Mass function for galaxies in the MegaZ-LRG (solid 
thick line) and the DR4 sample (dashed thin line). 

of the stellar mass function has apparently changed very 
little out to a redshift of 2 ~ 0.9. 

3.2 Radio Luminosity Function 

A standard technique for quantifying the rate of evolution of 
a population of galaxies is to compare the radio luminosity 
function of these objects at two different epochs. In this 
section, we will determine the evolution of the radio-loud 
AGN population by comparing the luminosity function of 
MegaZ-LRG radio galaxies with that derived from the SDSS 
DR4 sample. 

These functions were calculated for radio sources above 
a 1.4 GHz flux den sity limit of 3.5 mjy using the standard 
i/Vmax formalism (jSchmidtl Il968 !). We estimate Vmax for 
each galaxy as the maximum comoving volume within which 
it would be visible given the radio flux and optical magni- 
tude limits and the colour cuts used to define the MegaZ- 
LRG sample (see Section 2.1). We have included both K- 
corrections and e-volutionary corrections in this computa- 
tion. These corrections were described in the previous sec- 
tion. Radio luminosities are calculated using the formula 

loglo[Pl.4GHz] = logio[47rD^(z)Sl.4GHz(l + z)""'] (5) 

where Dl is the luminosity distance in the adopted cosmol- 
ogy, Si.4GHz is the measured radio flux, (l-|-z)"~^ is the stan- 
dard fc-correction used in radio astronomy, and a is the radio 
spectral index (5^ oc for which we ado pted a value of 

0.7, a s is usually assumed for radio galaxies (|Condon et al.l 
|2002| ). 

The normalization of the luminosity function requires 
knowledge of the precise intersection area of all three sur- 
veys: NVSS, FIRST and the SDSS DR4 photometric sample. 
We used the footprin t services of the Virtual Observatorj0 
l|Budavari et al.ll2006l ) to derive a total area of 6164.6 deg^. 
A further 8.3 deg'^ was subtracted to avoid noisy regions 
around bright NVSS sources, resulting in an effective area of 

^ http:/ /www. voservices.net/footprint 



6156.3 deg^. In the MegaZ sample, the maximum variation 
of the mean redshift as we go from our second luminosity bin 
to the highest luminosity one, is ~ 0.025. If we include the 
first bin, this difference increases to a value ~ 0.08. We nev- 
ertheless decide to keep the lowest luminosity bin to avoid 
introducing further luminosity cuts. 

In Figure 5 we plot the radio luminosity function mea- 
sured for the MegaZ-LRG catalogue and the results are also 
tabulated in Table 3. For comp arison we reproduce the lu- 
minosity function calculated bv I Sadler et al.l l|2007l ) for the 
2SLAQ LRG sample. There is very good agreement, but our 
new determination has considerably smaller error bars, be- 
cause our sample is very much larger. This is especially im- 
portant for determining the number density of radio sources 
with radio powers greater than logio(Pi.4GHz)= 26 W Hz"^. 
Also plotted in the figure are two determinations of the ra- 
dio luminosity function from SDSS DR4: the first without 
color restrictions, and the second restricted to red galaxies 
as defined by a luminosity dependent g — r color cut found 
by fitting bi-Gaussian functions to the distr ibution of q — r . 
The details of the method are explained bv lLi et al.l (|2006l ). 
For reference purp oses, we also show the RLF obtained by 
iBest et al.l l|2005al ) for the DR2 sample. 

The evolution of the luminosity function of a popula- 
tion of objects is often characterized in terms of simplified 
density or luminosity evolution models. In a pure density 
evolution model, the underlying luminosity distribution of 
the radio sources remains fixed, but the number density of 
the sources increases at higher redshifts. This might be the 
case if a larger fraction of galaxies harbour radio jets at 
higher redshifts, but the intrinsic properties of the jets do 
not evolve with time. In a pure luminosity evolution model, 
there is a fixed population of radio AGN with luminosities 
that dim as the radio sources age. 

In the two right-hand panels of Figure 5, we plot the 
ratios /pdens and /plum as a function of radio luminosity. 
/pDENS is the ratio of the comoving number densities of 
radio sources in the MegaZ and DR4 catalogues at fixed 
radio luminosity, while /plum is the factor by which the low 
redshift number densities should grow in order to match the 
values found at higher redshift. If either model is a correct 
description of the data, one should find constant / values as 
a function of radio luminosity. We find that neither model 
is a good description of the data. Both plots show a sharp 
increase near logio (Pi.4ghz)~ 25 W Hz"^. 

In Figure 5, we also plot the radio luminosity func- 
tion for steep-spectrum sources at z=0. 55 predicted by the 
luminosity /density evolution model of iDunlop fc Peacock 
((1990i ). This model is based on high frequency data, and has 
been shown to provide an accurate description of the radio 
source counts and redshift distributions. At the redshifts cor- 
responding to our MegaZ-LRG sample, the models were only 
strongly constrained above 10^® W Hz~^ due to the high flux 
limits of their radio samples, and were then extrapolated 
to the lower radio powers. We note that iDunlop fc Peacockl 
( 199d ') provide a range of different parametrizations in their 
paper, but all give very similar results at moderate redshifts. 
As ca n be seen, our data matches the iDunlop fc Peacock 
(|l990l ) model reasonably well. In partcular, the extrapola- 
tion to lower radio power appears to provide a reasonable fit 
to our data, but our estimate appears to be 0.2-0.3 dex lower 
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24 25 26 27 28 24.5 25.0 25.5 26.0 

log,o(L,.4GHz / W Hz-') log,o(L,,,GHz / * Hz"') 

Figure 5. The radio luminosity function at 1.4 GHz derived for MegaZ-LRG and DR4 samples (thick line with square and circle 
symbols, respectiv ely). An addi tional DR4 RLF restricted only to red sequence galaxies (triangles) i s also p lotted. Likewise shown are 
the calculations bv lSadler et al .' ('2007) for 2SLAQ (thin line and diamond symbols), and bv lBest et~aLl ll2005al) for SDSS DR2 (dashed line 
and cross es) . Results are in g o od ag reement at both redshift intervals. The luminosity function based on the luminosity /density evolution 
model of iDunlop fc Peacockl l|l99ol ) is also overplotted (dash-dotted line). The right panels show the factors by which luminosity and 
density should increase in order to match the DR4 curve at {z) ~ 0.14 with the MegaZ— LRG curve at (2:) ~ 0.55. No evolution scenarios 
are indicated by the horizontal dashed lines. Error bars are poissonian, except for the right upper panel where they are calculated via 
statistical bootstrapping of each sample. 



that that of iDunlop fc Peacockl l|l990l ) over the luminosity 
range from 10^^ to lO^'^ '^ W Hz'^ 

It is important to note that if radio galaxies fall out- 
side the MegaZ-LRG colour selection criteria, they will be 
missed in our radio luminosity function determination. Such 
an effect is visible for our low redshift luminosity function 
at radio powers above 10^® W Hz~^; the luminosity func- 
tion of bright sources in the red galaxy population is 0.1-0.2 
dex lower than for the DR4 sample as a whole. Neither the 
DR4 sample nor the MegaZ-LRG sample includes radio- 
loud quasars, which would likely increase the amplitude of 
the radio luminosity functions by ~ 0.1 dex above a radio 
luminosity of 10^^ W/Hz~^. Nevertheless, it is clear from 
our comparisons that a large fraction (i.e > 50%) of all ra- 
dio galaxies at 2 = 0.55 do reside in the LRG population. 
In the rest of the paper, we will be comparing radio AGN 
fractions at low and high redshifts, and we will show that 
such fractional estimates are much less sensitive to any in- 
completeness in our sample. 

3.3 Fraction of Radio-Loud AGN 

We are now in a position to investigate the relation between 
radio-loud AGN and their host galaxies. Over the redshift 
range covered by the MegaZ-LRG sample (0.4 < z < 0.8), 
radio sources have 1.4 GHz powers well above 10^'' W Hz'^ 



Table 3. The MegaZ-LRG luminosity function for radio— loud 
AGN at 1.4 GHz 



loglo(Pl.4GHz) 


logio(p) 


N 


(W Hz-1) 


(Mpc-3 (logioL)-i) 




24.45 


c- 10+0.02 
'^•^•^-0.02 


2218 


24.75 


5 Qn+O Ol 
^•■^^-0.01 


5438 


25.05 


K KO+O.Ol 
''•■^^-0.01 


3506 


25.35 


c OQ+0.02 


1742 


25.65 


-6 20+° °^ 
"•^'J-0.02 


840 


25.95 


6 54+0 03 
-D.04_Q Q3 


365 


26.25 


(■ 70+0.04 


185 


26.55 


7 14+O.O6 


87 


26.85 


7 p.o+0.07 
' •'^•-'-0.09 


41 


27.15 


7 7C-+0.12 


18 


27.45 


00+0.20 


3 


27.75 


_o 1O+0.24 
^•^•-'-0.60 


1 



(the faintest radio galaxies in our sample are 10^* * W 
Hz~^). It is well known that radio emission above this lu- 
minosity is rarely associated with star- formati on, but rather 
an AGN (jSadler et aLll2002l : lBest et al.ll2005al ). In addition, 
luminous red galaxies are photometrically defined so that 
the majority of star-forming galaxies are excluded, so we 
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Figure 6. The left upper panel shows the percentage of galaxies have radio luminosities above a given threshhold, as a function of 
stellar mass. Results are shown for the high redshift MegaZ-LRG sample (thick solid lines), and for the low redshift SDSS DR4 sample 
(thin dashed lines). The right upper panel shows the radio— loud fractions for the SDSS DR4 sample, but including a colour constraint 
of g — r ^ 0.8. Lower panels show the relative variation 5f between high redshift (MegaZ) and low redshift (DR4) radio-loud fractions 
as a function of stellar mass. Errors are poissonian and are propagated via standard error propagation theory. 



can safely assume that our sample is dominated by an AGN 
population. 

Figure 6 shows the fraction of MegaZ LRGs with ra- 
dio luminosities above diflterent threshholds, as a function 
of their stellar mass. The fraction of LRGs with radio AGN 
more luminous than lO^*'* W Hz~^ rises from ~0.04% for 
galaxies with stellar masses around 10^°-*Mq, to ~25% for 
the most massive LRGs with M. ~ IO^^'^Mq. Results for 
radio power cuts of lO^^ W Hz-\ lO^^'^ W Hz"! and 10^^ 
W Hz~^ are also included in the plot and a similar mass de- 
pendence is found. The dotted lines connecting light symbols 
on the plot show the corresponding dependence of the radio 
loud fraction on stellar mass for the low-redshift SDSS DR4 
sample. We find that there is an increase at higher redshifts 
in the fraction of of radio-loud AGN hosted by galaxies of 
all stellar masses. The increase is larger at higher radio lu- 
minosities and for galaxies with lower stellar masses. This 
is shown in more detail in the bottom three panels of Fig- 
ure 6, where we plot the relative ratio 5f between the high 
redshift and low redshift data {Sf = (/McgaZ — /dr,4)//dr4) 



as a function of mass and for our three radio luminosity 
thresholds. 

Once again it is important to recall that the DR4 
sample includes all galaxies regardless of colou r, whereas 
the Me gaZ sample is restricted to red galaxies. iBest et al] 
(l2005bl ) showed that the probability for a galaxy to host a 
radio-loud AGN was not influenced by its stellar population. 
In the top right panel of Figure 6 we compare the radio-loud 
AGN fraction versus mass relation for the full DR4 sample 
to the results obtained if the sample is restricted to galax- 
ies with g ~ r colours greater than 0.8. As can be seen, the 
results are almost identical. The effect of the colour cut is 
much smaller than the evolutionary trends we find between 
the low redshift and MegaZ LRG sample. 

To check whether the increase in the fraction of radio- 
loud AGN with mass is not simply the result of the fact that 
more massive galaxies harbour more massive black holes and 
hence produce more radio emission, we have also calculated 
the fraction of radio-loud AGN that are above a given limit 
in radio luminosity per unit stellar mass. For our sample. 
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Figure 7. The diflfcrcncc between the true radio-loud AGN fraction and the one calculated for the artifically rcdshifted samples is 
plotted as function of stellar mass. Results are shown for increasing luminosity cuts. The solid line shows the mean difference for the 25 
random samples, and the gray area is the 1-a scatter in this variation. Very few sources have masses higher than 10^^ '^ Mq. 



this is equivalent to choosing only those sources that pro- 
duce radio emission above a fixed fraction of the Edding- 
ton Umit. For reasonable values of the chosen limit (e.g. 
L/M* > 10^^'^WHz"^Mq^), the strong mass dependences 
persists, although the slope of the relation is somewhat shal- 
lower. 



3.4 Other systematic effects 

Photometric redshift codes have improved vastly over the 
years and the resulting redshift estimates arc useful for a 
broad range of studies. Even though the neural network pho- 
tometric redshift estimator used in the MegaZ catalogue has 
been throughfuUy calibrated using spectroscopic redshifts, it 
is important to test the senstivity of our results to photo- 
z errors. Such errors could translate into systematic offsets 
in our estimates of stellar mass. In addition, the errors in 
the Sloan u,g,r,i,z magnitudes increase with redshift, so 
one might worry that this would translate into apparent (as 



opposed to real) evolution of the radio loud fraction with 
redshift. 

In order to test for such effects, we artificially redshifted 
our SDSS DR4 sample from their observed redshifts near 
(z) ~ 0.14 out to the mean redshift of the MegaZ-LRG sam- 
ple, z = 0.55. To do this, we applied the kcorrect algorithm 
in conjunction with the e-corrections described above. We 
then added an error to the redshift using the distribution of 
Zphot — Zspec of 2SLAQ galaxies. In practice, for each MegaZ 
LRG, we extracted the 100 closest neighbors in g, r, i, z mag- 
nitude space in the 2SLAQ catalogue and we draw the red- 
shift error, Zphot Zspe-c-j at random from this set of galaxies. 
Magnitude errors were assigned to each DR4 object in the 
same way. Wo then ran our code for estimating stellar mass 
on this "artificially redshifted" data set, first applying just 
the redshift error, and second, incorporating errors on both 
the photometry and the redshifts. The whole procedure was 
repeated 25 times. 

The solid line in Figure 7 shows the difference between 
the mean fraction of radio-loud AGN measured in the 25 



Evolution of the Radio Loud Galaxy Population 11 




log,o(L,.4GHz / W Hz-^ log,o(L,.4GHz / W Hz-^ 

Figure 8. Integral bivariate radio luminosity-mass function in several bins of stellar mass (symbols) for high redshift MegaZ-LRG 
sources (thick solid lines) and low redshift SDSS DR4 sources (thin dashed lines). The bivariate function gives the fraction of radio— loud 
AGN brighter than a given radio luminosity, as a function of the luminosity. The right panel shows the ratios of the bivariate functions 
derived for the MegaZ-LRG sample to those derived for the SDSS DR,4 sample. These results show that radio sources evolve strongly 
across cosmic time depending on their stellar mass and radio luminosity. The dashed horizontal line marks the no— evolution scenario. 



artifically redshifted samples and the "true" value for the 
DR4 galaxies as a function of stellar mass. The difference 
5f is expressed in fractional terms, i.e. Sf = (Jkrtiflciai — 
Jtruc)/^truc. The shadcd contours respesent the la scatter 
among the 25 samples. It can be seen that small systematic 
shifts in the derived radio loud fractions do occur at the 
two ends of the stellar mass distributions. They also occur 
in our estimates of the fraction of the most luminous radio 
galaxies. On the whole, these effects are small (< 20%). Note 
also that 10^^'^ W Hz~^ sources are much rarer in the DR4 
catalogue than in the MegaZ sample, so the bottom right 
panel overestimates the size of the true effect at this radio 
luminosity and is likely to apply only to very much higher 
luminosity cuts at 2 = 0.55. 



3.5 Bivariate radio luminosity— mass function 

The first determination of the biva riate radio-optica l lumi - 
nosity function was carried out by lAuriemma et al.1 (|l977l ) 
for a sample of 145 galaxies o f type E and SO. and 
there were subsequent studies by I Sadler et al.l l| 19891 ) and 
iLedlow fc OwerJ (|l996l ).This (integrated) function measures 
the probability that a galaxy in a given optical absolute 
magnitude range hosts an AGN with radio power above a 
certain value. Best ct al. (2005b) calculated the bivariate ra- 
dio luminosity-stellar mass function of radio sources in the 
SDSS DR2, and parameterized the behaviour using a bro- 
ken power with characteristic luminosity P* = 2.5 x 10^*W 
Hz~^. A single fitting function was found to hold for for all 
masses below IO^'^Mq. 

Our MegaZ radio galaxy sample is o ver 100 times larger 
than the one of lAuriemma et "all (| 19771) and p icks up close 
to the break luminosity f ound by Best et al] l|2005bl ). Fol- 
lowing iBest et al.l (|2005bh . we will use stellar mass instead 



of optical luminosity. We can thus interpret our calculation 
as the fraction of radio-loud AGN of a given mass emitting 
above a certain certain radio power. Figure 8 shows the bi- 
variate functions at high and low redshift in mass bins from 
10" Mq to IO^^'^'^Mq. The most striking result shown in 
this plot is that for galaxies of all stellar masses, the high 
redshift fractions lie only slightly above the low redshift ones 
at faint radio luminosities, but there has been much stronger 
evolution at higher radio luminosities. This is shown in more 
detail in the right hand panel of Figure 8, where we plot the 
ratio of the radio-loud AGN fraction derived for the MegaZ 
LRGs to that derived for the DR4 galaxies, as a function 
of radio power. We show results for the four mass bins for 
which there is reasonable overlap between the two surveys 
(there many fewer galaxies with 12 < log M, < 12.25 in the 
DR4 sample, because the volume covered is much smaller). 
It is striking that for all the mass bins, the ratio remains 
approximately constant up to a radio luminosity of 10^^ W 
Hz~^ and then turns up sharply. As discussed in Section 1, 
lO'^^ W Hz~^ corresponds to the approximate dividing line 
between low excitation and high excitation sources, and also 
delineates the boundary between strongly-evolving and less 
strongly-evolving radio sources found in previous studies. 
The main new result shown in this plot is that the strongest 
evolution in the radio source population apperently occurs 
in the lowest mass galaxies. 

Another way of showing the same thing is to plot the 
integrated radio emission per unit stellar mass for galaxies 
of different stel l ar ma sses in the DR4 and MegaZ samples. 
iHeckman et al.l (|2004l ) used the volume-weighted distribu- 
tion of [OIII] 5007 line luminosity per unit black hole mass 
to study the average rate at which black holes are accret- 
ing and growing. They find that low mass black holes are 
currently growing at substantially higher rates than high 
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Figure 9. Integrated volume-weighted radio luminosity per unit 
stellar mass, as function of stellar mass, for MegaZ-LRG radio- 
loud AGN (solid line with square symbols) and DR4 sources 
(dashed line with diamond symbols). 



mass black holes. iBest et al.l l|2005bl ) extended this analysis 
by computing the integrated radio luminosity per unit black 
hole mass, demonstrating that radio and emission line lumi- 
nosities are prod uced in black holes of quite different masses 
(see Figure 11 in lBest et aLll2005bl) . 

In the present work, we do not have direct estimates 
of black hole mass, so we carry out the computation as a 
function of stellar mass. In Figure 9 we plot the volume- 
weighted radio luminosity per unit stellar mass, as a function 
of log Af* , for our high and low redshift samples. The plot 
is constructed by taking at each mass bin, the ratio of the 
integrated (and weighted) radio luminosity, to the integrated 
(and weighted) stellar mass in that bin. As can be seen, radio 
activity has been boosted at all masses at high redshift, but 
the boost factor is considerably higher for low mass galaxies. 



4 THE STELLAR MASSES AND COLOURS OF 
THE HOST GALAXIES OF RADIO-LOUD 
AGN 

It has long been known that powerful high redshift radio 
galaxies are hosted by very massive galaxies with predom- 
inantly old stellar populations. This conclusion arises from 
the fact that there is a well defined relation between the 
K-band magnitude of radio galaxies and redshift (the fa- 
mous K-z relation) that agrees with the predicions of a pas- 
sively evolving model galaxy of high mass (> lO^^A/0), and 
that there is relatively little scatter about this relation for 
the most luminous radio sources (e.g. iLiUv fc Longai^l 19851 : 
[jariiseFaDliool). 

So far we have phrased our analysis of the relationship 
between radio AGN and their host galaxies in terms of the 
bability for a galaxy of a given mass to host such an AGN. 
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Figure 10. Median stellar mass of radio-loud AGN in the 
MegaZ-LRG (solid line) and SDSS DR4 (dotted line) samples as 
function of radio power. Gray dots represent the actual values for 
each MegaZ galaxy and errors are calculated by the bootstraping 
technique. 



We have shown that this probability evolves most strongly 
in low mass galaxies. Note, however, that the fraction of 
radio loud AGN is a strongly increasing function of galaxy 
mass both ai z — 0.55 and at z = 0.14. This means that 
radio galaxies in massive elliptical hosts always dominate by 
number. 

This is shown in more detail in Figure 10 where we plot 
the median stellar masses of both the MegaZ and the SDSS 
DR4 radio AGN as a function radio luminosity. As can be 
seen, radio AGN are hosted by galaxies with median masses 
of ~ 3 X IO^'^Mq in both surveys. There is a weak increase 
in the median mass at higher radio luminosities. 

In Figure 11, we plot the colour distributions of the 
MegaZ radio AGN. One must of course bear in mind that the 
MegaZ sample is selected so as to occupy a restricted range 
in colour space, but we can nevertheless check whether there 
are any colour differences between different types of radio 
AGN. We split the sample into 3 different redshift bins and 
show results for low luminosity (< lO'^^'^ W Hz~^) and high 
luminosity ( > 10'^^'^ W Hz~^) radio so urces. As was found 
previously bv lSmith fc HeckmanI (|l989l ). there is a tendency 
for the more powerful radio galaxies to have slightly bluer 
colours at fixed redshift, but the effect is very weak. 

In summary, we conclude that the majority of both low 
luminosity and high luminosity radio AGN are hosted by 
massive (3 L*) galaxies with red colours. This is in good 
agreement with past work on this subject. 



5 SUMMARY AND DISCUSSION 

The main results of the present work can be summarized as 
follows; 

A catalogue of 14453 radio-loud AGN with 1.4 GHz 
fluxes above 3.5 mjy in the redshift range 0.4 < z < 0.8, 
has been constructed from the cross-correlation of NVSS 
and FIRST radio source catalogues with the MegaZ-LRG 
catalogue of luminous red galaxies. The vast majority of 
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Figure 11. Colour distribution of MegaZ-LRG in three differ- 
ent redshift bins, split into a high luminosity population of radio 
sources with logio (Pi.4GHz)^ 25.3 W Hz~^, and a low luminosity 
population with logio(Pi.4GHz)> 25.3 W Hz~^. 



the radio AGN are single component sources in both NVSS 
and FIRST. However, there is a signifiant fraction of objects 
without high S/N FIRST detections, which are required for 
accurate identification of the optical counterpart. We thus 
introduced a method for analyzing the FIRST radio maps 
around the candidate positions of these sources. This allowed 
us to dig deeper into the FIRST survey and use lower S/N 
detections to pinpoint the location of the host galaxy. 

We have presented a new determination of the luminos- 
ity function of radio-loud AGN at 2 0.55 that is in excel- 
lent agreement with other results in the literature, but with 
notably smaller error bars. By comparing our radio luminos- 
ity function at z ~ 0.55 to that derived for a large sample of 
nearby radio AGN from the SDSS DR4, we find compelling 
evidence for strong cosmic evolution of radio sources. The 
comoving number density of radio AGN with luminosities 
less thanlO^^ WHz-^ increases by a factor of ~ 1.5 between 
z = 0.14 and z = 0.55. At higher lumiosities, this factor in- 
creases sharply, reaching values ~ 10 at a radio luminosity 
of 10^'' W Hz-\ Neither a pure luminosity evolution nor a 
pure density evolution scenario provides a good description 
of the data. 

We then turn to an analysis of how the relation between 
radio AGN and their host galaxies evolves with redshift. The 
fraction of galaxies with radio luminosities above a given 
threshhold is a steeply increasing function of stellar mass at 
both z ~ 0.1 and at 2; ~ 0.55. The fraction of radio loud 
AGN increases with redshift and this increase is largest at 
the highest radio luminosities and also for lower mass galax- 
ies. We have also calculated the bivariate radio luminosity- 
mass function at z ~ 0.14 and z ~ 0.55. Its shape does not 
appear to depend on mass at either redshift, but there is 
strong evolution in the shape at the bright-end of the func- 
tion. The fraction of galaxies brighter than 10^* W hz"'^ de- 



clines with significantly shallower slope at higher redshifts. 
Within the range 10^^-^^ W hz^\ simple power law fits to 
the DR4 bivariate function in Figure 8 (frac oc 4ghz)' P'"''" 
duce slopes 7 of -1.84±0.21, -1.66±0.02 and -1.83±0.04 
for the first, second and third mass bins. In MegaZ, we ob- 
tain -1.05 ± 0.02, -1.16 ± 0.05 and -1.20 ± 0.02 for the 
same luminosity range and mass bins. 

In short, two main conclusions have emerged from our 
analysis: 

• There is a characteristic luminosity of ~ 10^^ W Hz~^ 
below which the radio source population appears to evolve 
only very weakly with redshift. Above this characteristic lu- 
minosity, there is strong evolution, with the most powerful 
radio sources undergoing the largest increase in co-moving 
number density. These results are in broad agreement with 
past studies. 

• The strongest evolution in the fraction of galaxies that 
host radio-loud AGN takes place in the lower mass galaxies 
in our sample. 

The most plausible explanation of these trends is that 
there are two classes of radio galaxy, likely associated with 
the high excitation/ low excitation "dichotomy") that have 
different fuelling/triggering mechanisms and h ence evolve 
in different ways. This has also been argued bv lTasse et all 
(2008), whose analysis of a sample 1% of the size of ours in 
the XMM-LSS region provided hints of similar evolutionary 
properties. As discussed before, it has been hypothesized 
that the class of low luminosity AGN is likely associated 
with massive (~ 10^^ Mq) dark matter halos with quasistatic 
hot gas atmospheres. The abundance of such halos is pre- 
dicte d to remain appro ximately constant out to redshifts 
~ 1 l|Mo fc Whitdl2002l ') and this is broadly consistent with 
the very weak evolution we see in this population. These 
"hot-halo" triggered sources at z ~ 0.5 would then produce 
radio emission that roughly compensates the radiative losses 
of the hot gas. 

The reason for the remarkably strong evolution in the 
comoving abundance of high-luminosity radio AGN is some- 
what more difficult to understand. As shown in Figure 8, 
this evolution is most dramatic in galaxies at the low stel- 
lar mass end of our sample. It is thus tempting to link this 
population with the strongly accreting (and evolving) pop- 
ulation of luminous quasars, whose space densities also in- 
crease by factors of more than 10 over the redshift range 
studied in this paper. Recent work has shown that quasars 
are hosted in dark matter halos of masses ~ 10^^ Mo, m- 
depende nt of redshift or the optical luminos ity of the sys- 
tem (e.g. lCroom et al.ll2005l : rShen et al.ll2007h . Only ~ 10% 
of optical quasars are radio-loud, however, and this raises 
the question as to whether there is simply a short radio- 
loud phase during the lifetime of every optical quasar, or 
whether optically luminous AGN and powerful radio galax- 
ies are triggered under different physical conditions. If, as 
in quasars, these sources are not regulated by cooling from 
a surrounding hot gaseous halo, but by other mechanisms 
like interactions/mergers, then we would not expect to find 
them in an equilibrum state like in low excitation AGNs. 

One way to shed further light on these matters is to 
compare the clustering properties of quasars to those of the 
luminous radio galaxies in our sample. If their clustering 
properties are identical, this would favour the hypothesis 
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that all quasars experience a radio-loud phase. Recently, 
iKauffmann et aU (2008) found that nearby radio- loud AGN 
were more strongly clustered than matched samples of radio- 
quiet AGN with the same black hole masses and extinction 
corrected [OIII] line luminosities. It will be interesting to see 
if the same conclusion holds at higher redshifts and for more 
powerful systems. This will be the subject of a future paper. 
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